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Abstract 

We apply HQET to semi-leptonic B meson decays into a variety of excited charm 
states. Using three realistic meson models with fermionic light degrees of freedom, 
we examine the extent that the sum of exclusive single charmed states account for 
the inclusive semi-leptonic B decay rate. The consistency of form factors with the 
Bjorken and Voloshin sum rules is also investigated. 



1 Introduction 



The role of semi-leptonic -B-decay in determining the CKM matrix element is well known. 
An important ingredient in this role is the understanding of the hadronic form factor, 
or Isgur-Wise (IW) function [|I|]. For B D**ei'e decays the heavy-light limit is an 
excellent starting point .Q Heavy quark effective theory (HQET) |^ is useful in organizing 
the descriptions of such processes. In this limit the wave functions and energies of the 
light degrees of freedom (LDF) are sufficient to define all IW functions [H. 

As the number of measured B decays increases, other questions related to semi-leptonic 
decay can be addressed. Among these are the B —>■ D* spectrum shape, and the ratio 
of D* to D rates. Of particular interest here are the inelastic processes B D**Xeh'e, 
where D** could be Di or D2 (or their radial excitations), and X are any non-charmed 
hadrons. Interest in these inelastic decays stems partly from experimental data |5|, 
for exclusive inelastic decays, and partly from the relation between the exclusive decay 
rates and measured and theoretical results for the inclusive semi-leptonic rate. 

An additional motivation for considering inelastic exclusive channels concerns sum 
rules. The Bjorken sum rule [0, H relates the slope of the elastic IW function (at the 
zero recoil point) to inelastic IW form factors describing S to P-wave semi-leptonic B 
decays. A less known sum rule involving the 5* to P-wave form factors has been derived 
by Voloshin 0. It is the analog of the "optical" sum rule for the dipole scattering of light 
in atomic physics. It is not immediately obvious whether a particular model yields form 
factors which are consistent with both of these two sum rules. 

Since the explicit IW functions are defined in terms of a relativistic LDF wave functions 
with a fermionic light quark, we have considered three appropriate hadronic models, based 
on the Dirac or the Salpeter equations. The models considered are the Dirac equation 
with scalar confinement, the Salpeter equation with time component vector confinement, 
and a relativistic fiux tube model. In each case the strong coupling constant and the 

^We use the symbol D** to denote any charmed meson. 
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various quark masses are chosen to best fit the heavy-hght data. The resulting form 
factor predictions are quite similar and the variety of different models used provides an 
assessment of the confidence one has in predictions of this sort. 

In Section ^ we briefiy review our previous exclusive form factor results, and in Section 
^ we outline the three Dirac-like hadronic models that are used for our numerical results. 
In Section ^ predictions for elastic and inelastic branching ratios into a single charmed 
hadron are compared with other results that can be found in the literature. A discussion 
of theoretical and experimental results of fractional semi-leptonic decay rates is found 
in Section ^ By considering fractional inclusive rates it is plausible that much of the 
uncertainty associated with the Vet value, quark masses, and QCD corrections cancels out. 
The consistency of form factors with the Bjorken and Voloshin sum rules is considered in 
Section ^. A summary of the results and our conclusions are given in Section |^. 



2 Isgur-Wise form factors and semi-leptonic B de- 
cays in the heavy- quark hmit 

In the ruQ — »• oo limit the angular momentum of the LDF decouples from the spin of 
the heavy quark, and both are separately conserved by the strong interaction. Therefore, 
total angular momentum j of the LDF is a good quantum number. For each j there are 
two degenerate heavy meson states (J = j ± |), and we can label states as J^. 



In HQET the covariant trace formalism |T^, [TT], [T2| is the most convenient for keeping 
track of the relevant Clebsch-Gordan coefficients and for counting the number of indepen- 
dent form factors. Using the notation of the lowest lying mesonic states are labeled 



as follows: C and C* denote Oi and li states (L = doublet), {E, E*) and (F, F*) 

2 2 

denote the two L = 1 doublets (Ot, it) and (it, 2t), respectively, and G and G* denote 

2 2 2 2 

ll and 23 states (L = 2 doublet). These states are described by 4 x 4 matrices, and 

2 2 

matrix elements of the heavy quark currents are calculated by taking the corresponding 
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traces |12|. In this way, one can write the decay rate for B D**eVe in the form 



-m 



duj 487r3 

Here uj = v ■ v' denotes velocity transfer, r** = mD**/TnB, and the function /** is given 
by 





(a;2-l)(l + rc)^ 


(2) 


fc*iuj,rc*) = 


{u + 1)[{lu + 1)(1 - rcO' + 4^(1 - 2iurc^ + r^.)] , 


(3) 


fEiuJ,rE) = 


(^2-l)(l-r^)2, 


(4) 


fE'{uJ,rE') = 


{u - l)[{uj - 1)(1 + rp^f + 4cu(l - 2u;rE* + r|0] , 


(5) 




^{00 -l){uj + - 1)(1 + rpY + uj{l - 2ujrp + rl)] , 


(6) 




|(cu -l){uj + l)^[{uj + 1)(1 - rp,y + 3cu(l - 2ujrp, + r|.)] , 


(7) 


fci^^rG) = 


^{uj - Ifiuj + l)[{iu + 1)(1 - rcf + toil - 2ujrG + r^)] , 


(8) 




'^{u - Ificu + l)[(c^ - 1)(1 + rc^Y + 307(1 - 2ujrG^ + r^,)] . 


(9) 



The above expressions can be found in |T^-[|T^. 

The only unknown quantity in the expression (|I|) is the appropriate IW form factor 
^**(ci;). Since these form factors cannot be calculated from first principles, one has to 
rely on some model of strong interactions. By comparing the wave function approach 
of |jT8| with the covariant trace formalism of [0, |Tl], |12|, and performing the necessary 



integrations in the modified Breit frame (as suggested in [T^), one finds the expressions for 
the unknown form factors in terms of the LDF rest frame wave functions and energies 0. 
These expressions include transitions from the ground state into radially excited states, 
and are given as 

2 

^c(^) = — —{jo{ar))oo 
UJ + 1 

2 

= / 2 1 0'iK))io 

y/uj^ — 1 



^g(^) = ^J^02(ar))2o 





C, C* transitions) , 


(10) 


(c- 


E,E* transitions) , 


(11) 


(c- 


F, F* transitions) , 


(12) 


(c- 


G, G* transitions) . 


(13) 



In the above formulae a is defined in terms of the initial {Eq) and final {E'^ energies of 
the LDF as 

a = (^. + ^;)\/^- (14) 
The expectation values appearing in (p!0|)- ([l3|) are defined as 



{FirT)-:- = [ r^rfr[/:,^:,(r)/^(r) +^:^.,(r)^^(r)]F(r) . 



(15) 

This follows from the form of the wave function in the Dirac-like models with spherical 
symmetry, 



-.k /^\^i-k, 



(16) 



where y^^^. are the usual spherical spinors, k = I {I = j + ^) or k = —I — 1 (/ = j — i), 
and a denotes all other quantum numbers. 



3 Modelling the light degrees of freedom 



In order to obtain a reasonable estimate of the model dependence of our results, we employ 
three qualitatively different realistic models to describe heavy-light mesons: the Dirac 
equation with scalar confinement (DESC), the Salpeter equation with vector confinement 
(SEVC), and the relativistic flux tube confinement (RFTC). All three models involve a 
short range Coulomb potential. With wave function of the form (|16|), it can be shown 



T9|, 20 1 that all three models satisfy a radial equation of the form 



{E^l -Ho- 



fHr) 



where 



m„ 



-rrin 



and 



k 1. 



:i7) 



(19) 



The 2x2 matrices L and I will be defined when the specific models are discussed. The 
numerical methods used to deal with these three models are described in [|19|, ^ . 



3.1 Dirac equation with scalar confinement (DESC) 

Scalar confinement is the only type of confinement potential that has correct sign of the 
spin-orbit coupling. In the Dirac equation it also yields linear Regge trajectories. This 
model also assumes a time component vector short range Coulomb interaction. Specifi- 
cally, we have 



1 



3r 



V 



3r 



(20) 



(21) 



The parameter values chosen to give an excellent fit (see Table |l|) to the heavy-light spin 



averaged data |2T| are 



rrir 



rrib 



a. 



O.SOOGey (fixed) , 
0A65GeV , 
1.357GeV , 
4:.693GeV , 
0.462 , 

0.284GeV^ (fixed) 



(22) 



The quality of fit is insensitive to the value of and the confinement tension b was 

chosen to yield the universal Regge slope |1T7||.FI 

^The slope of the Regge trajectories in the heavy-Ught case is expected to be exactly twice the slope 

, i.e. a'^^ = 2q;^^. The observed Regge slope for the light-light states is 



in the light-light case |19 



0.. 
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3.2 Salpeter equation with vector confinement (SEVC) 

The instantaneous version of the Bethe-Salpeter equation ^ (usually referred to as 
the Salpeter equation [|^) is widely used for the discussion of bound state problems. It 
is also equivalent [|23 to the so called "no-pair" equation 



which was introduced in 



order to avoid the problem of mixing of positive and negative energy states that occurred 
in the Dirac equation for the helium atom. A similar problem also occurs for a single 
fermionic particle moving in the confining Lorentz vector potential. For a very long time 
29| it has been known that there are no normalizable solutions to the Dirac equation in 



this case. 



It has been shown analytically for the heavy-light case |jl9[, and numerically for the 
case of fermion and antifermion with arbitrary mass ||3y, that in this type of model 
linear scalar confinement does not yield linear Regge trajectories. We have therefore used 
time component vector confinement with short range Coulomb interaction, even though 
it is well known that this model gives the wrong sign of the spin-orbit coupling. 

In this model matrices L and I are given by 



L 



2E2 



3r 







3r 



(23) 



(24) 



where 



and 



2Et 



mi 



D^D± . 



(25) 



(26) 



Again we fix the confinement tension b to yield the universal Regge slope, and choose our 
parameters as 



mu,d = 0.300GeV (fixed) 



m, = 0.598Ge1/ , 

= lA06GeV , 

rrib = A.lAlGeV , (27) 

as = 0.539 , 

b = O.U2GeV^ (fixed) , 

to obtain an excellent fit to the spin averaged heavy-light states, as shown in Table 0. 



3.3 Relativistic flux tube confinement (RFTC) 

In the RFTC model formalism for fermionic quark confinement is unusual in that the 
confinement is introduced into the kinetic rather than in the usual interaction term. The 
flux tube contributes to both energy and momentum, so it makes little sense to consider 
it as a "potential" type interaction. By a covariant substitution we add the tube to the 
quark momentum and energy. We may equivalently view this as a "minimal substitution" 
of a vector interaction field. The result nicely reduces to the Nambu string in the limit 
in which the quark moves ultra-relativistically. This physically motivated generalization 
of the potential model incorporates many aspects of QCD pO| . 

In this model the L matrix is the same as in (|23|) , while the interaction matrix I is 
given by 



/ 



(28) 



In the above T+ is defined as 



2 



[l-k) 



--Pt 



+ k 
--Pt 



(29) 



The flux tube energy and momentum, obtained by symmetrization of the classical expres- 
sions 



H, H, are defined by {{A, B} = AB + BA) 



±k 



Pf 



+ k 

a arcsm vj^ 
= a{r,F(^f)}, 



(30) 
(31) 



with (7^ = l/Jl-w2) 



1 ^arcsinf^ 1 , 



(32) 



Av± v± 

The only unknown operators in the above expressions are f J^. These are determined from 
the heavy-hght orbital angular momentum equation as in the spinless RFT model 



With the definition Wj. = Jp'^ + m?, these equations are pO 



= UWr^^ivij + air^Fivi)} 



jm V 



(33) 
(34) 



The numerical technique used to solve for v± is discussed in detail elsewhere p3| 
Theoretical predictions of the model with parameters 



rus 
rric 

h 

a. 



0.300 GeV (fixed) , 
0.580 GeV , 
1.350 GeV , 
4.685 GeV , 
0.181 GeV"^ (fixed) 
0.508 , 



(35) 



are shown in Table [1|. Again, the agreement with spin-averaged experimental masses is 
very good. 

We conclude this section by noting that all of the above models have been used for the 
calculation of the elastic IW form factor |1^, and the predicted IW functions {^ciyj)) 
were all quite consistent with the experimental data 



4 Branching ratios and comparison with other re- 
sults 

Our results for branching ratios obtained from the three different models discussed in 
Section ^ are shown in Table ||. We have assumed here that Vch = 0.040 and tb = 1.5ps. 
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Table ^ contains a comparison of our DESC calculation with calculation of Scora et al. 
(ISGW2 model) |3|, the one of Suzuki et al. (SISM model) |16|, with the QCD sum rule 
approach of Colangelo et al. (CNP model) [0, and with the model of Sutherland et al. 
(SHJL model) It is worth noting that results quoted for SISM and CNP are also 

obtained in the heavy quark limit. Table ^ contains ratios of partial widths for B decays 
into members of the same D** doublet. 

For calculation of branching ratios we have used experimental meson masses wherever 
possible. In those cases the only model dependent inputs were the appropriate IW form 
factors. For decays where the D** mass is unknown, we have used spin-averaged masses 
obtained in a specific model. Based on the available information on the splitting between 
D and D* (or Di and 02)-, one could estimate mass splitting in other excited doublets, 
and use that together with model dependent spin-averaged mass to obtain separate pre- 
diction for the mass of each member of that doublet. Meson masses obtained in this way 
could then be used in the calculation of the branching ratio for the corresponding decay. 
However, we have found that this procedure does not significantly affect the results. For 
example, using spin-averaged mass of 1974 MeV for D(1867) and Z}*(2009), instead of 
their experimental masses, in the case of DESC yields branching ratios of 2.235% and 
6.773% instead of 2.401% and 6.615%, which are given in Table H. For the higher states 
this effect is even less noticeable. 

Let us first discuss B ^ D and B ^ D* transitions. Recent results from CLEO p9| . 



BR{B- ^ D^e-Ue) = (1.95 ± 0.55)% , (36) 
BR{B- ^ D*^e-Ue) = (5.13 ± 0.84)% , (37) 

and results given in Tables ^ and ^ imply that all models we used, as well as the ISGW2 
and SHJL models, require Vcb slightly lower than 0.040. In our models values range from 
about 0.036 for DESC, to about 0.038 for RFTC and SEVC. On the other hand, ISGW2 
is consistent with 0.035, SHJL gives about 0.036, and SISM and CNP models agree with 
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Vcb of about 0.040. From Table |^ it can also be found that^ 



BR{B- 



D^e-jj, 



BR{B- D*^e-Ue) 



0.48 , for ISGW2 

0.34 , for SISM 

0.33 , for CNP 

0.31 , for SHJL 



(3^ 



On the other hand, our calculation with three different models yields (see Table |) 



BR{B- 



BR{B- 



0.35 , for RFTC 
0.36 , for DESC 



(39) 



0.35 , for SEVC 

The results quoted in ( |36D and (^) imply an experimental ratio of 

BR(B- D^e-Ve) , . 
) ^ = 0.38 ±0.17. 40 

BR{B- D*^e-Ue) 

It is interesting to note that ratio of polarization states of D and D* is 0.33. 

Individual contributions of P-wave j = | states to the total semi-leptonic decay rate 
is another interesting point. From Table § it can be seen that the total semi-leptonic 
branching ratio for B ^ Di and B ^ D2 is expected to be 

2 



BR{B- ^ Dle-Ue)+BR{B 



0.84 




0.040 


0.69 




0.040 


0.79 




0.040 



l.SOps 



l.SOps ■ 



l.SOps 



0.040 



for RFTC 
for DESC • (41) 
for SEVC 

and signif- 



cb 



0.040 



l.SOps 



These results are slightly larger than the ISGW2 result of 0.65 
icantly disagree with SISM, CNP, and SHJL models (0.20, 0.37, and 0.46 
respectively). However, as one can see from Table ^ our ratios of these two P wave 
decays are qualitatively different from the one obtained in ISGW2 and SHJL, and agree 



■^It is worth noting that the calculation of used form factor definitions which are not consistent 
with the covariant trace formalism. 
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with SISM and CNP models. We find 



BR{B- D\e~Ve) 
while the other models yield 



0.60 
0.63 
0.60 



for RFTC 
for DESC 
for SEVC 



(42) 



BR{B- 



Dfe-u, 



2.00 
0.70 
0.50 
1.54 



for ISGW2 
for SISM 
for CNP 
for SHJL 



(43) 



Again, it is interesting to note that the ratio of number of polarization states of Di and 
D2 is 0.6. It remains to be seen whether this discrepancy between our results (which are 
obtained in the heavy-light limit) and ISGW2 and SHJLQ models can be explained with 
large ^ effects [|0[. 

As already mentioned, within the HQET framework the only model dependent input 
for the decays C C,C* and C ^ F, F* are the unknown IW functions. For these decays 
the uncertainty introduced by using a particular model is about 10%. For other decays, 
D** mass is not known, so that in the calculation of decay rates we used predictions of 
a particular model. Therefore, one should expect larger discrepancies between different 
models. From Table |^ one can see that this is indeed the case. However, results obtained 
from the three different models are not significantly different. For example, for the decays 
C E,E* uncertainty the introduced by using a specific model is about 20%. Also, as 
one can see from Table ^, the ratios of the two exclusive decay widths for members of the 
same doublet are all consistent, which is the consequence of the application of HQET. 



■*In |p8| one can also find results obtained in the heavy quark limit. These are in general much smaller 
than our results, but the ratios of partial widths for the B decays into the members of the same D** 
doublet agree much better with our predictions. 
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5 Fractional semi-leptonic decay rates 

The exclusive decay rates discussed earlier suffer from a variety of theoretical oversim- 
plifications. Some of the things which were not taken into account are QCD corrections, 
spectator effects, and deviations from exact heavy quark symmetry. In addition, there 
are several parameters which need to be specified before definite predictions can be made. 
Among these are the CKM parameter Vcb, the 6-quark lifetime, and the quark masses. 

Many of the above problems can be reduced by considering fractions of the inclusive 
b —>■ ceUe rate. In particular, Vcb exactly cancels. Also, since the sum of the exclusive 
rates equals the inclusive rate, and since the inclusive calculation is structurally similar to 
the exclusive ones, there should be some cancellation of the QCD, spectator corrections, 
and heavy quark mass dependence. Since the inclusive rate has been measured, one can 
directly compare these fractional predictions with experiment in several cases. 



The inclusive spectator model decay rate for b — > ceue is [41, 42 



where 43 



I{x, 0, 0) = (1 - -8x + x^) - 12a;2 \nx . (45) 

If for the moment we ignore the oversimplifications of the above inclusive model and 
assume Vcb = 0.040, the 6-quark lifetime of Tf, = 1.5 x lO^^^s, and the quark masses of the 
three realistic models given in (P^D, (P^), and (P5|), we find a total branching ratio 

10.27 \Vcb/OMOf {tb/1.50ps)% , for RFTC 

10.31 |Kf,/0.040|^ {tb/1.50ps)% , for DESC • (46) 

10.53 |K=fe/0.040|^ {tb/1.50ps)% , for SEVC 
The experimental branching ratio, 



BR{b — > ceue 



BR{B- Xe'Ve) = (10.49 ± 0.46)% , (47) 

is in excellent agreement with the above numbers. However, one should keep in mind that 
the predicted value is very sensitive to the choice of Vcb, U, and quark masses. 

13 



One plausibly assumes that ratio of the exclusive branching ratios to the inclusive one 

BR{B D**eue) 



R* 



(48) 



BR{b cez7e) ' 

will be more accurate than either of these separately. We first apply this idea to B ^ D 
and B ^ D* decays. From (^), ( P7| ) and ( ^Tj) we find experimental fractions, 



0.19 ±0.06 



0.49 ±0.10 . 



(49) 
(50) 



From Table ^ one can see that our three models predict 



R 



0.20 , for RFTC 
0.23 , for DESC 
0.20 , for SEVC 



(51) 



and 



r)th 



0.58 , for RFTC 

0.64 , for DESC • (52) 
0.58 , for SEVC 

The predicted values are reasonably consistent with measurement in all cases. 

The fraction of semi-leptonic decay into final states other than D or D* is by (^9]) and 



m 



^expp** ^^Yier than D and D* 



1 _ + 



0.32 ±0.16 . 



(53) 



From Table ^ we see that three models discussed in this paper imply 



K^{D** other than D and D* 



1-{R'I^ + R%) 
0.22 , for RFTC 
0.13 , for DESC 
0.22 , for SEVC 



(54) 
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It is interesting to observe that single excited charmed states alone are nearly consistent 
with accounting for the entire inclusive semi-leptonic decay fraction. As a more direct 
way of seeing this note the total fractional percentage in Table 0. The predicted fraction 
into all D** states lies at 89% or above for the three models discussed in this paper. 



6 Consistency with sum rules 
6.1 The Bjorken sum rule 

The Bjorken sum rule Q relates the derivative of the elastic form factor to the values 
of inelastic S- to P-wave form factors at the zero recoil point. In our notation, 



(55) 



Since the S- to P-wave form factor normalizations at the zero recoil point are not fixed 
(as is the case of the elastic form factor ^c)i but instead depend on the energies and wave 



functions of the LDF |T^, it is not immediately obvious that form factors obtained 
from the three different models will also satisfy the Bjorken sum rule. In particular, we 
observe from (|55D the manifestly valid inequality 

1 



(56) 



On the other hand, it follows from (p!0[) that 15 



(57) 



Combined together, these two bounds imply 



1 



(i)/ 



2 1 
> - 

- 4 



(5^ 



One may ask here whether the above two bounds for —^^(1) are consistent, or can one 
devise a model for which the S- to P-wave form factors in (^) come to less than ^ 0. 

With this in mind we have evaluated numerically both sides of (|55|) in order to check 
self consistency of the three models used in this paper. In the sum of the right-hand side 
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of ([55|) we included the lowest P-waves plus the first two radial excitations. For the Dirac 
equation with scalar confinement, with parameters given in (p2D , the direct evaluation 
yields 

-ec(l)- 0.86, (59) 
while the sum rule approach (the right-hand side of (p^D) gives 



^^(1) = 0.25 + (0.535 + 0.026 + 0.004 + ...)- 0.82 . (60) 



For the Salpeter equation with vector confinement, with parameters given in (|27|), we 
obtained after direct evaluation 

-ec(l)-l-26, (61) 

and the sum rule result was 

-e^(l)^ 1.03. (62) 
Similarly, the relativistic fiux tube model confinement, with parameters given in (|35D , 



yields after the direct evaluation 

-^^1)^1.14, (63) 

while the sum rule approach gives 

-ec(l)-l-09. (64) 

It is also worthwhile noting here the CLEO result for derivative of the renormalized form 
factor ^ 

-^"^(1) = 0.84 ±0.12 ±0.08 . (65) 

Since all the terms on the right-hand side of (^) are positive definite, and since we 
are neglecting nonresonant contributions to final states containing a pion (even if those 
were small), one might argue that in a self-consistent model result for —,^^(1) obtained 
by direct calculation should be smaller than the one obtained in the sum rule approach. 
Indeed, this is what happens in all three of the heavy-light models used in this paper. 
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Furthermore, for all three models (and especially for the DESC and RFTC) the sum rule 
is very close to being saturated by the resonant contributions. 

However, we must also mention at this point that a similar calculation was also per- 
formed in p3]. These authors find (also in the heavy quark limit) 



-ec(l) = l-28, (66) 

after the direct calculation, and 

-^^(1) = 0.25 + 0.21 + ... ^ 0.46 , (67) 

in the sum rule approach (using only the lowest P-wave mesons) . The above result shows 
that the sum rule is far from being saturated by resonances. The difference between the 
two approaches was in explained as being mainly due to nonresonant contributions 
to final states containing a pion. 



6.2 The Voloshin sum rule 

The Voloshin sum rule ||^ is the analog of the "optical" sum rule for the dipole scattering 
of light in atomic physics. In terms of our form factors and energies of the LDF [Ej:,** = 
rriD** — rric), it can be written in the form 

I = J E(%^ - 1) eg^l) + ^ E(%^ - 1) #(1) - A . (68) 

Here i^jj** and i?j;^L denote energies of the LDF of the j = \ and J = | P-wave mesons, 
respectively. In this sum rule the dependence on the unknown charm mass (or the LDF 
energy) is even stronger than in the case of the Bjorken sum rule, where it is contained 
only implicitly through form factors. Therefore, one could expect that model calculations 
used with (|68|) are less satisfactory than in the case of (|55|). Still, one could again make 
the same arguments as in the case of the Bjorken sum rule, and conclude that any self- 
consistent model calculation of the right-hand side of (|68D should yield result smaller than 
0.5. 
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To evaluate the right-hand side of (|68|), we have used the spin-averaged energies of 
the LDF in D and D* and other D** doublets (because the spin-averaged experimental 
masses were used in numerical calculations which determined quark masses). As before, 
we have used the lowest P-waves plus the first two radial excitations for both, j = \ and 
j = I mesons. 

For the Dirac equation with scalar confinement (DESC), with parameters given in 



(P^), we find 

A = 0.395 + 0.033 + 0.007 + . . . ~ 0.44 , (69) 

which is smaller than the predicted value of 0.5. On the other hand, the Salpeter equation 
with vector confinement (SEVC), with parameters given in (P7D, results in 

A ~ 0.63 , (70) 



while the relativistic flux tube model confinement (RFTC), with parameters given in (|35|) , 
yields 

A ^ 0.60 . (71) 

These two results are about 20% higher than the predicted value of 0.5. Again, we empha- 
size that in the case of the Voloshin sum rule numerical results are much more dependent 
on the unknown charm quark mass and other parameters of the model. Unfortunately, 
we are not aware of any other model which has used the Voloshin sum rule as a test of 
self-consistency, so that we cannot compare our result with the literature. Nevertheless, 
we do have to say that two out of three heavy-light models used here (SEVC and RFTC) 
appear to be inconsistent as far as Voloshin sum rule is concerned. Based on the above 
calculations, one might also argue that the DESC model predictions are the most reliable 
of all the results presented in this paper. 
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7 Conclusion 



We have examined the role of semi-leptonic B decay into higher charmed mesons. Within 
a HQET framework we have evaluated branching ratios for B —>■ D**eue, where the D** 
are all 5*- and P— wave mesons, wave mesons with j = |, and some of their radial 
excitations. Our numerical calculations are based upon three realistic models. In each case 
a light fermion interacts with a fixed source. A short distance vector Coulomb interaction 
is used, and at large distances the fermion is confined by a scalar, time-component vector, 
or a flux tube.0 

An important aim here has been to determine how much of the semi-leptonic decay rate 
ends up as one of the higher charmed resonances. Considered as a fraction of the inclusive 
b ceUe rate we find that between 90 and 95% of B mesons decay semi-leptonically into 
a single excited charmed resonance. 

Another area which we have explored in this paper is the consistency of form factors 

obtained from three different models with the Bjorken and Voloshin sum rules. One 

might suspect that there may be an inconsistency between the Bjorken expression for 

the slope of the Isgur-Wise function at the zero recoil point, and the expression for the 

same quantity which is obtained in the wave function approach and depends only on the 

elastic {S- to S-wave) transition. The two methods yield different upper bounds for the 

IW slope. However, using numerical values for the lowest two S- to P-wave form factors 

(and also for the derivative of elastic form factor) at the zero recoil point, we find that 

all three models yield results consistent with the Bjorken sum rule. The Voloshin sum 

rule provides another sensitive test of model calculations. There, the dependence on the 

unknown charm quark mass is even stronger than in the case of the Bjorken sum rule, 

and one might expect that ^ effects play an even more significant role. All three of the 

heavy-light models used in this paper give results which are slightly higher than the sum 
^Results based on spinless models such as semi-relativistic quark model (also used in (l^), or spinless 
relativistic flux tube model (used in [^) are not significantly different from the results based on the 
models discussed in this paper. 
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rule prediction. Nevertheless, given the large qualitative differences between models itself, 
and the fact that their predictions are very similar in all cases, we believe that all the 
results presented in this paper are trustworthy. 

We have also compared our results with other calculations available in the literature. 
We find significant disagreements with and in ratios of decay widths for B decays 
into the members of the same D** doublet. On the other hand, our results are in general 



significantly larger than the ones obtained in 16 and [37 
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TABLES 



Table 1: Heavy-light states. Spin- averaged experimental masses are calculated in the 
usual way, by taking | (|) of the triplet and j (|) of the singlet mass for the S{P) waves). 
Theoretical results are obtained from the three heavy-light models discussed in the text. 
Model parameters are given in (H) (DESC), (|7D (SEVC), and (H) (RFTC). Theoretical 
errors with respect to the spin-averaged experimental masses are shown in parentheses. 



State 






h 

rC 


Mass 
(MeV) 


RFTC 

(^iViev ) 


DESC 
(^iViev ) 


SEVC 
(^iViev ) 


cu, cd quarks 
















L>(1867) 


C 


07 
















2 


-1 


IS (1974) 


1981(+7) 


1977(+3) 


1980(+6) 


D*(2009) 


C* 


i 1 

2 














Di(2425) 


F 


It 
















2 


-2 




2439(-7) 


2444(-2) 


2440(-6) 


L>^(2459) 


p* 


2t 
2 




'P2 










cs quarks 


















D,(1969) 


C 


07 




'So 












2 


-1 


IS (2076) 


2071(-5) 


2074(-2) 


2072(-4) 


D;(2112) 


C* 


17 

2 














D,A2535) F 


1 + 
J- 3 






IP (2559) 












2 


-2 


2564(+5) 


2560(+l) 


2564(+5) 


0*2(2^73) F* 


2t 
2 














bu, bd quarks 
















5(5279) 


c 


07 




'So 


1^ (5314) 










2 


-1 


5316(+2) 


5313(-1) 


5316(+2) 


5*(5325) 


c* 


17 

2 




's. 










bs quarks 


















fis(5374) 


c 


07 




'So 












2 


-1 


1^ (5409) 


5407(-2) 


5410(+1) 


5407(-2) 


fi:(5421) 


c* 


17 

2 




's. 
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Table 2: Exclusive partial widths for decays B D**eve obtained from the three different 

2 



models discussed in the paper. F is given in units of 

2 ^ 



0.040 



BR is in units 



of 



cb 



0.040 



l.SOps • 



while the ratio R 



BR{B^D**ev^) 
BR(b^cePe) 



is given in [%]. Numerical values of 



BR{b — > cez/g) for a particular model can be found in 



State 




RFTC 






DESC 








SEVC 








r 


BR 


R 


r 


BR 


R 




r 


BR 


R 


c 


07 

2 


9.026 


2 


057 


20.03 


10.54 


2.401 


23.29 


9.420 


2.147 


20.40 


c* 


17 

2 


26.09 


5 


946 


57.89 


29.03 


6.615 


64.14 


26.80 


6.108 


58.03 


E 


0+ 

2 


0.211 





048 


0.468 


0.303 


0.069 


0.670 





204 


0.047 


0.442 


E* 


It 

2 


0.299 





068 


0.663 


0.419 


0.096 


0.926 





289 


0.066 


0.626 


F 


It 

2 


1.383 





315 


3.069 


1.161 


0.265 


2.565 


1 


294 


0.295 


2.802 


p* 


2t 

2 


2.307 





526 


5.119 


1.854 


0.423 


4.097 


2 


152 


0.490 


4.659 


G 


17 

2 


0.016 





004 


0.036 


0.023 


0.005 


0.051 





016 


0.004 


0.035 


G* 


27 

2 


0.016 





004 


0.036 


0.023 


0.005 


0.051 





016 


0.004 


0.035 


G2 


07 

2 


0.225 





051 


0.499 


0.067 


0.015 


0.148 





215 


0.049 


0.466 


c; 


17 

2 


0.460 





105 


1.021 


0.131 


0.030 


0.289 





437 


0.100 


0.946 


E2 


ot 

2 


0.009 





002 


0.020 


0.018 


0.004 


0.040 





022 


0.005 


0.048 


E* 


It 

2 


0.011 





003 


0.024 


0.025 


0.006 


0.055 





028 


0.006 


0.061 


F2 


It 

2 


0.068 





016 


0.151 


0.029 


0.007 


0.064 





085 


0.019 


0.184 


F* 


2t 

2 


0.101 





023 


0.224 


0.045 


0.010 


0.099 





128 


0.029 


0.277 


total 


40.22 


9.17 


89.24 


43.67 


9.95 


96.49 


41.11 


9.37 


89.01 
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Table 3: Our results with the DESC model for B D**ei>e compared to predictions of 



ISGW2 13 , SISM |T|, CNP |3, and SHJL ||. T is given in units of 

2 ^ „ 



0.040 



and BR in units of 



Vet 



0.040 



l.SOps ■ 



btat 


e 


ISGW2 


SISM 


CNP 


SHJL 


DESC 


U 


tP 


r 


BR 


r 


BR 


r 


BR 


r 


BR 


r 


BR 




Ui 

2 


12.53 


2.860 


7.800 


1.778 


7.616 


1.736 


9.478 


2.160 


10.54 


2.401 


G 


i 1 

2 


26.12 


5.950 


23.17 


5.279 


23.39 


5.331 


30.646 


6.980 


29.03 


6.615 


771 

h, 


Ui 

2 


0.316 


0.072 


0.118 


0.027 


0.272 


0.062 


0.295 


0.067 


0.303 


0.069 




i 1 

2 


0.316 


0.072 


0.154 


0.035 


0.381 


0.087 


0.421 


0.096 


0.419 


0.096 


F 


2 


1.896 


0.432 


0.363 


0.083 


0.544 


0.124 


1.232 


0.281 


1.161 


0.265 


p* 


2t 

2 


0.948 


0.216 


0.517 


0.118 


1.088 


0.248 


0.800 


0.182 


1.854 


0.423 


G 


I3 

2 


not 


given 


0.001 


0.000 


not 


given 


0.022 


0.005 


0.023 


0.005 


G* 


23 

2 


not 


given 


0.001 


0.000 


not 


given 


0.003 


0.001 


0.023 


0.005 


C2 


07 

2 


0.000 


0.000 


0.071 


0.016 


not 


given 


0.579 


0.132 


0.067 


0.015 


CI 


17 

2 


0.632 


0.144 


0.172 


0.039 


not 


given 


not { 


jiven 


0.131 


0.030 


E2 


ot 

2 


not 


given 


not 1 


given 


not 


given 


not { 


jiven 


0.018 


0.004 


E* 


It 

2 


not 


given 


not 1 


given 


not 


given 


not { 


pven 


0.025 


0.006 


F2 


It 

2 


not 


given 


not 1 


given 


not 


given 


not J 


pven 


0.029 


0.007 


F* 


2t 

2 


not 


given 


not J 


given 


not 


given 


not { 


pven 


0.045 


0.010 


total 


42.76 


9.74 


32.36 


7.38 


33.29 


7.59 


43.49 


9.91 


43.67 


9.95 
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Table 4: Ratios of partial widths for the B decays into the members of the same D** 
doublet obtained from the three different models discussed in this paper (second to fourth 
columns), compared to ISGW2 |§, SISM |T§, CNP [0, and SHJL |38[ results. 



Doublet 


RFTC 


DESC 


SEVC 


ISGW2 


SISM 


CNP 


SHJL 


C/C* 


0.35 


0.36 


0.35 


0.48 


0.34 


0.33 


0.31 


E/E* 


0.71 


0.72 


0.71 


1.00 


0.76 


0.71 


0.70 


F/F* 


0.60 


0.63 


0.60 


2.00 


0.70 


0.50 


1.54 


G/G* 


1.00 


1.00 


1.00 


not given 


0.98 


not given 


7.00 


G2/G2 


0.49 


0.50 


0.49 


0.00 


0.41 


not given 


not given 


E2/E2 


0.82 


0.72 


0.79 


not given 


not given 


not given 


not given 


F2/F* 


0.67 


0.64 


0.66 


not given 


not given 


not given 


not given 
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